Evaluation of entropy driven jet symmetry transitions

Here we determine whether entropy drives planar turbulent jets into round turbulent jets. Determining when a jet flow transitions from one symmetry to the next is an important but incompletely resolved problem. The constructal view argues that the transition between symmetries of jet flows is governed by the minimization or maximization of entropy. Here we explore whether entropy increases with the transition of a planar turbulent jet into a round turbulent jet and whether entropy maximization (or minimization) predicts the same location of symmetry transition as velocity matching. We find that entropy considerations presented do not predict this transition.

Here we determine whether entropy drives planar turbulent jets into round turbulent jets. Determining when a jet flow transitions from one symmetry to the next is an important but incompletely resolved problem. The constructal view argues that the transition between symmetries of jet flows is governed by the minimization or maximization of entropy. Here we explore whether entropy increases with the transition of a planar turbulent jet into a round turbulent jet and whether entropy maximization (or minimization) predicts the same location of symmetry transition as velocity matching. We find that entropy considerations presented do not predict this transition.
Why turbulent jets change symmetry is an open question. For example, rectangular turbulent jets become turbulent circular jets a suitable distance away from the orifice. That transition is particularly interesting because it is not induced by interaction with a surface (e.g., circular turbulent jets transforming into radial wall jets). However, why this transition occurs in the absence of a surface remains unclear.
Bejan et al. 1 argue that observation of turbulent planar jets morphing into round turbulent jets is an example of evolution. They support this argument by contending that this change in jet symmetry increases jet mixing as an example of "nature's arrow", phenomena traditionally associated with entropy. Yet, Bejan et al., are careful to not use the term entropy, arguing from a constructal view that entropy must only remain positive for irreversible processes, and use velocity matching alone to determine the location of the symmetry transition without reference to entropy development. Heitor Reis argues that minimization or maximization of entropy are direct corollaries of the constructal view, suggesting that entropy (or a rate thereof) may play an essential role in the transition between jet symmetries 2 .
Inspired by this potential explanation of jet symmetry transition, here we explore whether entropy increases with the transition of a planar turbulent jet into a round turbulent jet and whether entropy maximization (or minimization) predicts the same location of symmetry transition as velocity matching. We find that entropy considerations presented do not predict this transition.

Model
Consider, as do Bejan et al. 1 fully immersed turbulent jets in single component fluids (e.g., water) of infinite domain. The jet comprises the system of interest. Entropy is related to internal energy, pressure, and volume through where T is the absolute temperature, S is the specific entropy, U is the specific internal energy, p is pressure, and V is the specific volume (inversely proportional to the density ρ). Application of the total derivative (D/Dt) indicates within the control volume defined by the jet, where t is time. Internal energy conservation may be expressed as where q is the heat flux, τ is the deviatoric stress tensor, and v is the local fluid velocity vector 3 . Substitution and simplification for isothermal jets in the absence of mass sources or sinks (e.g., suction) 4 then requires with μ here specified as the local turbulent viscosity (i.e., for moderate to low viscosity fluids) and Γ is the magnitude of the rate of strain tensor 3 . The square indicates that viscous dissipation is always positive as expected. We now evaluate the double dot product for the planar and round free turbulent jets.
Consider then a planar turbulent jet issuing in the stream wise z direction from a slot of length L and width W, where L ≫ W, with spreading in the x direction in Cartesian coordinates. With velocity and derivatives in the y direction (perpendicular to z) vanishing (i.e., v y = d/dy = 0), simplified at constant density with continuity as where the subscript indicates a component in the direction indicated. The error in this expression relative to the three-dimensional version is O((W/L) 2 ) where O is the order operator, so that when L ≫ W the error introduced by the simplification is quantitatively small. Because the velocity is turbulent, it may be divided into both timesmoothed averages and fluctuations relative to the average with v i = < v i > + u i , with < u i > = 0, where i = {x, z} and < … > indicates temporal averaging. Applying the averaging operator as defined by Deen 3 returns which shows that in addition to two terms for time-averaged velocities, two additional terms with instantaneous velocities also contribute. The time-average planar jet velocity profile has long been known to be self-similar sufficiently far away from the issuing orifice, with the velocity in the z direction as where the second equality defines v z,max for a planar jet, h p is the velocity-decay constant for the planar jet, U o is the average nozzle velocity, W is the nozzle width (not half width), and f p is the velocity distribution for the planar jet as a function of both x and z 5 . The impact of uncertainty from the self-similarity assumption is evaluated a posteriori. To first order approximation, f p is a Gaussian distribution that depends on a similarity parameter η p = x/β p z, where β p represents the jet spread of the planar jet. Mass conservation through the continuity equation specifies with use of symmetry to set the lower limit of integration. Rajaratnam 5 gives the turbulent viscosity using the Prandtl-Boussinesq eddy viscosity formulation as where k p is a constant. The fluctuating velocities may be approximated to first order as u i ~ (u i 2 ) 1/2 . Experimental data reported by Heskestad 6 and Rajaratnam 5 Fig. 1a, I p f = 3.65, which remains nearly an order of magnitude smaller, suggesting that instantaneous terms remain secondary quantitatively.
In contrast, consider a round turbulent jet issuing in the z direction from a nozzle of diameter d o with r radially transverse to the jet's axis of symmetry in cylindrical coordinates with velocity and derivatives in the azimuthal direction vanishing (i.e., v θ = 0, d/dθ = 0, (∇v) θθ = 0), Because the velocity is turbulent it may be divided into both time-smoothed averages and fluctuations relative to the average with v i = < v i > + u i , with < u i > = 0, where i = {r, z} and < … > indicates temporal averaging 3 . Applying the averaging operator as defined by Deen 3 returns (13)

Results and discussion
The purpose of this article is to determine whether entropy governs the transition between hydrodynamic jet symmetries. The remainder of this article first determines whether the transition is related to maximization or minimization of entropy (mutually exclusive options) and second determines whether predicted transitions are in quantitative agreement with experiment. Figure 2 presents viscous dissipation for both jets as a function of distance from the nozzle. Because viscous dissipation is positive for both jets, entropy is generated with both jets as anticipated from the second law of thermodynamics for irreversible processes. The figure shows that viscous dissipation decreases with distance from the nozzle. The round jet starts at a higher viscous dissipation rate but then decays more rapidly than the planar jet that starts at a lower viscous dissipation rate. If maximizing entropy generation (i.e., maximizing momentum mixing) were to have been the driver of the transition in jet symmetry then one would expect the viscous dissipation of the planar jet to exceed that of the round jet at short distances but the round jet to exceed the planar jet at large distances. However, the figure shows exactly the opposite. Therefore, the transition is not accompanied by an increase in entropy generation; entropy maximization does not govern the transition. This finding leaves the question of whether entropy minimization quantitatively predicts the transition as discussed below.
The difference in exponent on distance from the nozzle indicates an intersection of the curves. www.nature.com/scientificreports/ jet, 3β p zL (neglecting spreading in y), versus the radial spread of a round jet, π(3β r z) 2 /4, suggests a transition from planar to circular jets but at a location different than that proposed by Bejan, et al. 1 at This coefficient takes on values of 5.1-6.1 using values cited here. Experimentally, Deo et al. 7 and Sfeir 8 typically do not identify a simple transition between jet symmetries but find departure from Eq. (9) at meaningfully smaller values than adoption of Eq. (20), between which neither power law holds. Deo et al. 9 find these two transitions at z transition /L = 1.1 and 2.2, respectively. Re-analysis of the results of Sfeir 8 finds the first transition at z transition /L = 1.8-3.4 and the second transition at z transition /L = 5.9-8.0. When the two transitions coincide, they occur at z transition /L = 7.9-9.5. In either case, entropy considerations appear to not do a better job of predicting the transition between jet symmetries than the simple velocity matching of Sfeir and Bejan, et al.
The analysis above assumed the jets to be self-similar. Sangras and Faeth 10 show that the uncertainty associated with use of the self-similar assumption is approximately less than 30% at 10 nozzle diameters for a round turbulent nonbuoyant jet and negligible by 30 nozzle diameters. Quinn's experiments with slot jets show deviation from similar streamwise velocity profiles only at or below 5 equivalent diameters 11 . For conditions in Fig. 2, the transitions occur at > 30 equivalent nozzle diameters. Even with modest uncertainty (note the log scale in Fig. 2a), the conclusion that entropy considerations do not predict jet symmetry transitions does not change.

Conclusion
In summary, determining when a jet flow transitions from one symmetry to the next is an important but incompletely resolved problem. We find that entropy considerations presented do not predict the transition from planar to round jet. Therefore, the assertions of Sfeir 8 remain in force. Velocity matching "should be used only as a rough indication of the location of transition", and the "mechanism of transition" from planar to cylindrical symmetry "is not due to a simple redistribution of momentum" 8 . More advanced approaches that evaluate the stability of jet structures may be helpful, including nonlinear, weakly nonlinear, or Lyapunov-type stability approaches.